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, Abstract 
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I If neutrino masses and mixings are suitable to explain the atmospheric and 

solar neutrino fluxes, this amounts to contributions to FCNC processes, in par- 
ticular fj, — > e, 7. If the theory is supersymmetric and the origin of the masses is 
a see-saw mechanism, we show that the prediction for BR(/i — ► e, 7) is in gen- 
p ■ eral larger than the experimental upper bound, especially if the largest Yukawa 

coupling is 0(1) and the solar data are explained by a large angle MSW ef- 
fect, which recent analyses suggest as the preferred scenario. Our analysis is 
bottom-up and completely general, i.e. it is based just on observable low-energy 
data. The work generalizes previous results of the literature, identifying the 
dominant contributions. Application of the results to scenarios with approxi- 
mate top-neutrino unification, like SO(10) models, rules out most of them unless 
the leptonic Yukawa matrices satisfy very precise requirements. Other possible 
ways-out, like gauge mediated SUSY breaking, are also discussed. 
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1 Introduction 



In the pure Standard Model, flavour is exactly conserved in the leptonic sector since one 
can always choose a basis in which the (charged) lepton Yukawa matrix, Y e , and gauge 
interactions are flavour- diagonal. If neutrinos are massive and mixed, as suggested by 
the observation of atmospheric and solar fluxes 0, this is no longer true and there exists 
a source of lepton flavour violation (LFV), in analogy with the Kobayashi-Maskawa 
mechanism in the quark sector. Unfortunately, due to the smallness of the neutrinos 
masses, the predicted branching ratios for these processes are so tiny that they are 
completely unobservable, namely BR(/i — > ej) < 1CT 50 0. 

In a supersymmetric (SUSY) framework the situation is completely different. Be- 
sides the previous mechanism, supersymmetry provides new direct sources of flavour 
violation in the leptonic sector, namely the possible presence of off-diagonal soft terms 
in the slepton mass matrices ( m e H ) .; an d trilinear couplings A\- ||. There 

are in the literature very strong bounds on these matrix elements coming from requir- 
ing branching ratios for LFV processes to be below the experimental rates [[|. The 
strongest bounds come precisely from BR(/i — > ej). 

The most drastic way to avoid these dangerous off-diagonal terms is to impose 
perfect universality of the (mj^)--, (ml ) , AS, matrices, i.e. to take them proportional 
to the unit matrix. Then, of course, they maintain this form in any basis, in particular 
in the basis where all the gauge interactions and Y e are flavour- diagonal, so that we are 
left just with the previous tiny Kobayashi-Maskawa LFV effects. In principle, there 
is no theoretical reason to impose the universality constraint, and that is the so-called 
supersymmetric flavour problem. However, there are some particular supersymmetric 
scenarios, such as minimal supergravity, dilaton-dominated SUSY breaking or gauge- 
mediated SUSY breaking, that give rise to the desired universal soft terms (at least 
with much accuracy). In any case, the universality assumption is obviously the most 
conservative supposition that one can make about the form of the soft breaking terms 
when analyzing flavour violation effects; and this will be our starting point. 

It turns out, however, that even under this extremely conservative assumption, if 
neutrinos are massive, radiative corrections may generate off-diagonal soft terms. The 
reason is that the flavour changing operators giving rise to the (non-diagonal) neutrino 
matrices will normally contribute to the renormalization group equations (RGEs) of 
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the (m 2 )^-, (my., and Afj matrices, inducing off-diagonal entries. 

The most interesting example of this occurs when neutrino masses are produced 
by a see-saw mechanism 0. Then, above a certain scale M, there are new degrees 
of freedom, the right-handed neutrinos, which have conventional Yukawa couplings 
with a matricial Yukawa coupling Y u . Normally, Y e and Y u cannot be diagonalized 
simultaneously. Working in the flavour basis where Y e is diagonal, the off-diagonal 
entries of Y u (more precisely, the off-diagonal entries of Y+Y„) drive off-diagonal 
entries in the previous (rrv^)^, an d matrices through the RG running ||. 

Actually, the supersymmetric see-saw is a extremely interesting scenario for a num- 
ber of reasons. First, the see-saw mechanism is probably the most convincing and 
economical mechanism to naturally produce tiny neutrino masses (certainly it is the 
most popular one). Second, a ncm-supersymmetric see-saw suffers from a serious hi- 
erarchy problem. Perhaps this point has not been sufficiently acknowledged in the 
literature and we would like to emphasize it. It is common lore that the Standard 
Model (SM) presents a generic gauge hierarchy problem. Namely, the Higgs mass 
acquires quadratically-divergent radiative corrections, which, under the usual interpre- 
tation of the renormalization process, are naturally of the size of the scale at which 
new relevant physics enters. In a see-saw scenario the problem becomes more acute. 
The presence of very massive fermions (the right-handed neutrinos) coupled to the 
Higgs field produces logarithmically-divergent radiative corrections to the Higgs mass, 
which are proportional to M 2 . This can be immediately noticed from the usual ex- 
pression of the 1-loop effective potential, where the neutrino mass eigenstates which 
are essentially right-handed neutrinos contribute as g^M 2 !^ 2 ./? 2 (log ^ — |) + ... (/z 
is the renormalization scale) M . In the supersymmetric framework, this problem 
is automatically cured by the presence of the right-handed sneutrinos with a similar 
mass. Consequently, it is fair to say that supersymmetry is the most natural arena to 
implement the see-saw mechanism. 

This is precisely the scenario that we will consider in this paper, computing BR(/i — > 
e^) in the most general case compatible with the observational (atmospheric and solar) 
indications about neutrino masses and mixing angles. This subject has been addressed 
in previous works, especially in refs. (other related works can be found in 

ref.[0). Let us briefly comment on the differences between the treatment of these 
works and the present one. The authors of ref. |§ adopt a "top-down" point of view, 
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i.e. they start with some sets of particular textures for Y e and Y v , coming from more 
or less well motivated symmetries, e.g. family symmetries, which can produce sensi- 
ble low-energy neutrino mass matrices, Ai u , at the end of the day. Then, they study 
the corresponding predictions for BR(yU — > erj). On the other hand, the authors of 
ref . ||,|l(],|ll]] adopt a "bottom-up" point of view, i.e. they start with the experimental 
data about M. u at low energy, finding out a form of Y v able to reproduce these data, 
and then computing the corresponding BR(/i — > &y) . Both approaches are interesting, 
but not general, as they are related to particular choices of Y v . As a consequence, 
their results are not conclusive about what can we expect for BR(/i — ► ej) based on 
actual neutrino data. In our case, we also follow a "bottom-up" approach, i.e. we start 
with the existing experimental information about A4 U . Then we find the most general 
textures for Y v and Ai (the Majorana mass matrix of the right-handed neutrinos) con- 
sistent with that information. In our opinion, this is by itself an interesting issue. The 
textures turn out to depend on some free parameters, but very few (as we shall see, the 
textures used in ref.|ltj correspond to a particular choice for these parameters). Next, 
we compute BR(/i — > ej) and compare it with present and forthcoming experimental 
bounds. The approach can be considered as a generalization of that by Hisano et al. 
in ref. [||,[nj. As a matter of fact, we identify the dominant contributions, which were 
not considered in ref. PJlOU as a consequence of the particular texture used there. 

The paper is organized as follows. In sect. 2 we determine the most general form 
of Y v and Y+Y^, compatible with all the phenomenological requirements. In several 
subsections we specialize the general formulas previously obtained for relevant physical 
scenarios. In sect. 3 we review the way in which the RG-running induces off-diagonal 
terms, and the relevance of the latter for — > lj,^y processes. In sect. 4 we apply the 
information of sects. 2, 3 to compute the predictions on BR(/j — > lj, 7), with particular 
focus on BR(yU — > e, 7). These predictions can be understood analytically using some 
approximations, for the sake of clarity in the discussion. Again, in several subsections 
we present and discuss the results for different physical scenarios. In sect. 5 we offer 
a summary of the predictions on BR(/.; — > obtained in the previous sections. In 
sect. 6 we present our concluding remarks. Finally, we include an Appendix with the 
most important formulas that we have used for the full computations, specially RGEs 
and expressions for BR(/j — > 
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2 General textures reproducing experimental data 



The supersymmetric version of the see-saw mechanism has a superpotential 

W = Wo - \v c rMv c r + u c R T Y u L ■ H 2 , (1) 

where Wq is the observable superpotential, except for neutrino masses, of the preferred 
version of the supersymmetric SM, e.g. the MSSM. The extra terms involve three 
additional neutrino chiral fields (one per generation; indices are suppressed) not charged 
under the SM group: Um (i = e, /i, r). Y u is the matrix of neutrino Yukawa couplings, 
Li (i = e,fj,,r) are the left-handed lepton doublets and H 2 is the hypercharge +1/2 
Higgs doublet. The Dirac mass matrix is given by hid = Y U (H%). Finally, M. is a 3 x 3 
Majorana mass matrix which does not break the SM gauge symmetry. It is natural to 
assume that the overall scale of A4, which we will denote by M, is much larger than 
the electroweak scale or any soft mass. Below M the theory is governed by an effective 
superpotential 

W eff = W + ^(Y U L ■ H 2 ) T M-\Y V L ■ H 2 ), (2) 

obtained by integrating out the heavy neutrino fields in ([]]). The corresponding effective 
Lagrangian contains a mass term for the left-handed neutrinos: 

5C = ~v T M u v + h.c, (3) 

with 

M v = m u T M ^ = Y/M-'Y^) 2 , (4) 

suppressed with respect to the typical fermion masses by the inverse power of the large 
scale M. It is convenient to extract the Higgs VEV by defining the k matrix as 

« = M v /(Hl) 2 = Y v T M- l Y u , (5) 

where (H®) 2 = v 2 = v 2 sin 2 f3 and v = 174 GeV. The experimental data about neutrino 
masses and mixings are referred to the M. v matrix, or equivalently k, evaluated at low 
energy (electroweak scale). In this sense, it should be noted that eqs.(f|, |^) are not 
defined at low energy but at the "Majorana scale", M. Therefore, in order to compare 
to the experiment one has still to run k down to low energy through the corresponding 
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RGE0. Alternatively, following the bottom- up spirit, one can start with the physical 
k matrix at low energy and then run it upwards to the M scale, where eqs.(^, [5]) are 
defined. This has been the procedure we have followed in all the numerical calculations 
to be presented in the next sections. 

Working in the flavour basis in which the charged-lepton Yukawa matrix, Y e , and 



gauge interactions are flavour-diagonal, the k matrix is diagonalized by the MNS |T3 
matrix U according to 

U t kU = diag(Ki, k 2 , k 3 ) = D K , (6) 
where U is a unitary matrix that relates flavour to mass eigenstates 

(7) 

It is possible, and sometimes convenient, to choose Ki > 0. Then, U can be written as 

U = V ■ diag(e-^ /2 , e-^l\ 1) , (8) 

where <fi and 0' are CP violating phases (if different from or tt) and V has the ordinary 
form of a CKM matrix 
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(9) 



On the other hand, one can always choose to work in a basis of right neutrinos where 
M. is diagonal 

M = diag(.Mi, M 2 , M 3 ) = D M , (10) 



with M.i > 0. Then, from eqs.(|, |]) 



D K = U T Y/D M -iY v U = U'Y^D^jD^^YM 



T : 



where, in an obvious notation, D^j = +\/Da- Multiplying both members of the 
eq.(|TTJ) by D^—rr from the left and from the right, we get 



(12) 



2 For hierarchical neutrinos, with which we will be mainly concerned, this running does not sub- 
stantially affect the texture of k, and thus the mixing angles, ft does modify, however, the overall size 
of k, which has to be taken into account. 
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whose solution is D ^tj^tY V U D ^-rj = R, with R any orthogonal matrix (R can be 
complex provided R T R = 1). Hence, in order to reproduce the physical, low-energy, 
parameters, i.e. the light neutrino masses (contained in D K ) and mixing angles and 
CP phases (contained in U), the most general Y v matrix is given by 

Y„ = DvM RD^U + (13) 

So, besides the physical and measurable low-energy parameters contained in D K and 
U, Y v depends on the three (unknown) positive mass eigenvalues of the righthanded 
neutrinos and on the three (unknown) complex parameters defining R. We will see, 
however, that in practical cases the number of relevant free parameters becomes dras- 
tically reduced. Let us stress that eq . (p~3|) , which is central for us, has to be understood 



at the M scale and in the above-defined basis, i.e. the flavour basis in which the 
charged- lepton Yukawa matrix, Y e , the right-handed Majorana mass matrix, Ai, and 
the gauge interactions are flavour-diagonal. The extension of eq.(^) to other choices 
of basis is straightforward, taking into account the transformation properties of the 
Y e , Y u , M. matrices under changes of basis. For example, if one works in a basis 
of vr where the Ai matrix is non-diagonal, then the most general Y v matrix reads 
Y u = UljD ^jjRD , where Um is such that Ai = U* M R> jv^m- Let us also mention 
that eq.flTBD is valid as well for the non-supersymmetric see-saw mechanism. 

As commented in sect. 1, in ref.[|H| a particular choice for Y v was taken, which 
corresponds to take R — 1 in eq.([l3|). This is equivalent to assume that there exists 
a basis of Lj and vm in which Y v and Ai are simultaneously diagonal (though not 
Y e ). To see this notice that if R — 1, then Y u , as given by eq . flllf ) , can be made 
diagonal by rotating Lj with the £7-matrix. This hypothesis may be consistent with 
certain models, namely when all the leptonic flavour violation can be attributed to the 
sector of charged leptons, but clearly it is not general, and it is inconsistent with other 
scenarios. 

Another special choice of R occurs when the Y v matrix given by eq. ( |13"D has the form 
Yj, = WDy, where Dy is a diagonal matrix and W is a unitary matrix. From the see- 
saw formula (|5[), W is related to D M and D Y through D^j^t = W*(D y -ikD y -i)W + . 
Then, by rotating with the V^-matrix, Y e and Y v get simultaneously diagonal, 
while Ai gets non-diagonal, more precisely M. = DyK~ x Dy. In other words, all the 
leptonic flavour violation can be attributed to the sector of right-handed neutrinos, in 
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contrast to the previous situation. It is important to remark that both are very special 
situations, not generic at all. 

On the other hand, as we will see in the next section, the rate for U —>■ processes 
just depends on the Y+Y,, matrix (to be precise, on the matrix element (Y+Y„)jj) 
rather than on Y u , since this is the quantity that enters the RGEs of the off-diagonal 
soft parameters. From eq. (|i~3"D Y+Y„ has the general form 

Y+Y, = UD^- K R + D M RD^U + (14) 

Notice that Y+Y,,, and therefore eq.([L4]), do not depend on the i/R-basis, and thus on 
the fact that Ai is diagonal or not. 

In order to illustrate the use of eqs. ( |13| , 
that often appear in the literature 

2.1 z/^'s and vrs completely hierarchical 

In this case Ki < /t 2 C and M.\ <C M. 2 <C so that D K and D M can be 

approximated by 

D K ~ diag(0,K 2 ,K 3 ), D M ~ diag(0,0,^ 3 ) (15) 

where {v^k 2 ) 2 ^ (AmJ) s „i, {vl^f ^ (Am^) ffitm . Then, from eq.([TJ), 

(Y u ) tJ ~ sfM^R^Ufj (16) 

So, for a certain set of low-energy physical quantities (given by ki and U), Y u just de- 
pends on two (generally complex) parameters, and ^/AA 3R33. Let us remark 
that, although the {Y u )ij and (Y u ) 2 j entries are proportional to \/ M.\ and V ' M. 2 re- 
spectively, and thus suppressed, they are normally relevant if one wishes to reconstruct 
the neutrino mass matrix, k, from Y„, Dm through the see-saw equation (||). Still, they 
are not relevant for the Y+Y u entries, and thus for the predictions on BR(/j — > 
which is our main goal. 

It is interesting to notice that this scenario also includes the possibility of an "inverse 
hierarchy", i.e. ~ diag(A^3, 0, 0), or any other ordering. The reason is that one can 
permute the entries of Dm by redefining the z/r's through an appropriate rotation. In 
particular, one can pass from an inverse hierarchy to an ordinary one through a rotation 
in the 1-3 "plane": u R = R p v' R with R p = {(0,0,1), (0,1,0), (-1,0,0)}. Then, 
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14), let us consider some interesting scenarios 



Y v -»■ YJ = D' VM RlRD v JJ + with D' M = R T p D M R p . Hence, the general formula 
and its particularization ( |T6| ) remain identical with the modification R — > R' = R p R. 
Since R runs over all possible orthogonal matrices, so R' does. So, both starting points 
(inverse or ordinary hierarchy) lead to the same complete results. 

In many models one can impose an extra condition coming from the unification of 
Yukawa couplings. In particular, in simple 5*0(10) models, usually one of the Y+Y u 
eigenvalues (the largest one) coincides with the top Yukawa coupling \Y t \ 2 at the uni- 
fication scale, Mx (in other more elaborate models, it coincides with 3|Yt| 2 ) [[L4[]. We 



will denote \Yq\ this maximum eigenvalue. Now, from eq.(|16|), the largest eigenvalue 



of Y+Yy is |Y | 2 ~ M3 (|i?32| 2 ^2 + |-R33| 2 ^3), which (after running up to Mx) would 



32 < 



be identified with |Y t (Mx)| . This eliminates one real parameter, leaving just R33/R, 
plus one phase, e.g. phase(i?33), as the independent parameters. Actually, as men- 
tioned above, the U — ► processes just depends on the Y+Y„ matrix, which does 
not depend on the previous phase. Therefore, one is just left with one (complex) 
parameter, R^/R?^- 

A convenient way to visualize this is to take the following parametrizationP] of R 

' C2C3 — C1S3 — S1S2&3 S1S3 — C1S2C3 
R = c 2 s 3 C1C3 - S1S2S3 -S1C3 - C1S2S3 I , (17) 

V $2 S1C2 C1C2 



where #1, 82, 83 are arbitrary complex angles. Then, from eq. (|16|) 
Y„ 





Y 


c 2 




2 K 2 + 




^3 


C 2 | 



/ 

IL- 



which depends just on the angle §1 and the phase #r; Y+Y v depends just on 8\. The 
previous equation should be understood at the M. scale; the prefactor \Y \ represents 
the square root of the maximal eigenvalue of Y+Y^ at the .M-scale. 

In the typical case in which everything is real, very used in the literature, the 
&i angles are real, and hence Y u and Y+Yj, depend on a single real parameter, 8\. 



3 The parametrization ([l7]) does not count all the possible forms of R, as it does not include 
"reflections". However, for the scenario considered in this subsection it is general enough, since the 
uncounted forms of R do not produce different textures of Y„ . 
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Therefore, in spite of the additional degrees of freedom introduced by the lack of 
knowledge about R, this scenario is quite predictive. 

In the previous discussion, it has been assumed that at least one of the two entries, 
R32 or i?33, is different from zero. However, it may happen that R32 = R33 = 0. 



Then, K4, Aii, M.2 cannot be neglected anymore in the general equation eq.([13|). It is 
worth-noticing that going to the basis where the Dm matrix is 'inverse hierarchical', i.e. 
vr — ► RpVR, Dm — > diag(A^3, -M-i, Mi), this possibility translates into R12 = R13 = 0, 
which includes the particularly simple case R = 1. 

2.2 z/x's hierarchical and z/r's degenerate 
In this case, D K and Dm have the form 

D K — diag(0, k,2, 1^3), D M ^dmg{M,M,M) (19) 

where, again, (w|k 2 ) 2 ~ (Am 2 ) soi , (f^) 2 — (Am 2 ) atm . Then, from eq.(^3|), 

(Y„) - ~ VMRuVZPt (20) 

So, in principle, for a certain set of low-energy physical quantities (given by K q and U), 
Y v depends on one real parameter (Ai) and three complex ones defining the R matrix. 
Notice, however, that in the real case, where R + = R T , Y+Y u takes the simple form 

Y+Y, = MUD K U + , (21) 

which just contains one free parameter, Ai. This can be fixed by a unification condition, 
similarly to the previous case. Here |Y"o| 2 — — |^t| 2 (where the identification must 
be understood at the Mx scale). Therefore, this is also a very predictive scenario. 
Notice that here Y+Y^ does not depend on the form of R, so in this case the choice 
R = 1 is indeed general enough to study k — ► lj, 7 processes. 

2.3 z/^'s quasi- degenerate 

In this case D K ~ diag(/«i, K3), with K\ ~ k 2 ~ k% = k. Then, it is logical to 
assume that M. has degenerate eigenvalues, otherwise a big conspiracy would be needed 
between Y v and M. so that K\ ~ k 2 ~ ft 3 in eq.(g). Hence 

D K ~ diag( ), D M ~&a,g{M,M,M) . (22) 
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Consequently, eq.(|20|) and eq.(^l|) (for R real) hold in this case too. Since D K ~ k1, 
from eq.(^T|) we notice 



Y+Y„ = MUD K U + ~ |F | 2 1 , (23) 

where |Y" | 2 = A4k 3 is the largest eigenvalue of the Y+Yj, matrix. This means that for 
quasi-degenerate neutrinos, with R real, we expect small off-diagonal entries in Y+Y;, 
and thus suppressed k — > lj,^y rates. We will be more quantitative about the size of 
those entries in subsect. 4.3. 

The scenario of partially degenerate neutrinos will also be addressed in that section. 

3 RG-induced lepton flavour violating soft terms 

and li —> lj , 7 

In addition to the supersymmetric Lagrangian, which contains the gauge interactions 
plus the part derived from the superpotential (Q), one has to consider the soft breaking 
terms (gaugino and scalar masses, and trilinear and bilinear scalar terms) coming from 
the (unknown) supersymmetry breaking mechanism. In a self-explanatory notation, 
they have the form 

-Acft= \ m2 L) i M L i + ( m e a ) ii g » e flJ + ( A eij ejuHiLj ■,+ h.c.) + etc., (24) 

where we have written explicitly just the soft breaking terms in the leptonic sector, 
namely scalar masses and trilinear scalar terms. All the fields in the previous equation 
denote just the corresponding scalar components. As explained in sect. 1, concerning 
flavour violation the most conservative starting point for £ so ft is the assumption of 
universality, which corresponds to take 

{ m i) ij = m o 1 > { m2 e R ) i:j = m o A eij = A o Y eii , (25) 

so that working in the Li and em basis where Y e is diagonal, the soft terms do not 
contain off-diagonal (lepton flavour violating) entries. 

As mentioned in sect. 1, this is only strictly true at the scale where universality 
is imposed, e.g. M x in GUT models. Below that scale, the RGEs of the soft terms, 
which contain non-diagonal contributions proportional to Y+Y^, induce off-diagonal 
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soft termsQ These contributions are decoupled at the characteristic scale of the right- 
handed neutrinos, M. As noted e.g. in ref. [H], and as it is clear from eqs. (pTHTTl) in 



the Appendix, this mechanism is not efficient to generate off-diagonal entries in the 
( m ea) ma t r i x ) since the corresponding RGE does not contain ~ Y+Y^ terms. More 
precisely, in the leading-log approximation^, the off-diagonal soft terms at low-energy 
are given by 

04, - ^(3mg + ^)(Y+Y,),. log ^ , 
«)<, " ' 



{A e ) l3 ~ AoY k (Y+Y„)« log ^ , (26) 

where i ^ j and is the Yukawa coupling of the charged lepton l t . 
The amplitude for a — > lj, 7 process has the general form 

T = e%m u i(j af) qP (A L P L + A R P R ) h , (27) 

where q is the momentum of the photon, Prl = (l±7s)/2 and Al {Ar) is the coefficient 
of the amplitude when the incoming U lepton is left (right), and thus the lj lepton is 
right (left). The expressions for Al, Ar can be found in the literature (see Appendix 
for more details). The corresponding branching ratio is given by 

BR{l t - i il7 ) = ^ (\A L \ 2 + \A R \ 2 ) (28) 

Notice that the U —>■ ^,7 process violates the chiral symmetries Ir^ — ► e %a HR^. Under 
the assumption of initial universality, this symmetry is violated just by the Yukawa 
coupling of the corresponding lepton. Thus, Al {Ar) is proportional to the mass of 
the righthanded lepton involved, i.e. lj R (/^). Since m 2 j <C mi 2 , then \A L \ 2 <C \A R \ 2 . 
This can also be checked by inspection of the diagrams that contribute to U — > 
taking into account eq . (|26|) . It is interesting to note that the angular distribution of lj 
with respect to the polarization direction of the incoming lepton (PjJ depends on the 
relative magnitudes of the coefficients Al, Ar. The relation \Al\ 2 <C \Ar\ 2 predicted 

4 We will pay attention just to the contribution that appears in the MSSM, i.e. the one oc Y+Yj,. 
This is a minimal possibility. In more complicated models, there can be additional sources of flavour 
violation in the RGEs, leading to larger predictions for BR(/^ — > 

5 We use the leading-log approximation through the text in order to make the results easily under- 
standable. Nevertheless, the numerical results, to be exposed below, have been obtained by integrating 
the full set of RGEs. 
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Figure 1: Feynman diagrams contributing to U — > lj, 7. vx and lx with X = 1, • • • , 3 (4, • • • , 6) 
represent the mass eigenstates of "left" ("right") sneutrinos and charged sleptons, respectively. x~ A , 
A = 1,2, denote the charginos, whereas x^, A = 1, • • • , 4, denote the neutralinos. 




Figure 2: Dominant Feynman diagrams contributing to U — > ^,7 in the mass-insertion approxima- 
tion. Li are the slepton doublets in the basis where the gauge interactions and the charged-lepton 
Yukawa couplings are flavour-diagonal, \A denote the charginos and neutralinos, as in Fig. 1. 



by the see-saw mechanism, and which does not hold in general in other models of 
generation of neutrino masses, gives a characteristic (1 — P^cosO) distribution that 
could be measured by future experiments [|1^ . 

We have taken into account all contributions to BR(7j — > Zj, 7), using the general 
expressions given in the literature, in particular from ref. |J, as explained in the Ap- 
pendix. These have not been obtained by using the mass- insert ion approximation, 
but by diagonalizing all the mass matrices involved in the task, i.e. those of (left 
and right) sleptons, charginos and neutralinos. The diagrams have the form shown 
in Fig. 1. The precise form of BR(/j — > that we have used in our computations 
is a rather cumbersome expression, given in the Appendix. However, for the sake of 
the physical discussion it is interesting to think in the mass-insertion approximation 



to identify the dominant contributions. As discussed in ref. |]10| , these correspond to 
the mass-insertion diagrams enhanced by tan/3 factors. All of them are proportional 
to m 2 Lij , and have the generic form shown in Fig. 2. So, in all cases 
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DR(/, /,. - i ^\A R \ 2 ~ ^^%^tan^ , ,29) 



where m| represents supersymmetric leptonic scalar masses. 

From eq.(JZUD it is clear that BR(/j — > depends crucially on the quantity raf^ -, 



which at low energy is given by the integration of its RGE, as explained above. In the 
scenario examined in this paper (i.e. supersymmetric see-saw), and in the leading-log 
approximation, ra 2 Li - oc (Y+Y„)^-, as given by eq.(|26|), so 



a 



■■■•> 



2 



BR(k -> l j: 7) ~ ^— | ^(3ml + A%) log ^ (Y+Y„)<, 8 tan 2 (30) 



The Y^Y U matrix, which is therefore the crucial quantity, has been evaluated in sect. 2 
in the general and in particular cases. Next, we apply those results to the computation 
of BR(Z f 7 ), 

4 Experimental data on neutrinos and predictions 

for BR(/, ^,7) 

We study here the expectable predictions for Zj — > Zj,7 from the existing experimental 
information on neutrinos. Let us mention that the present and forthcoming experimen- 
tal bounds on U — ► make \i — *■ e, 7 the process with higher potential to constrain 
theories involving flavour violation. The present limit is BR(ii — > e, 7) < 1.2 x 1CT 11 
[W\ and it is going to improve in the near future (year 2003) until ~ 10~ 14 [fT7|| . Hence, 
we will pay special attention to /J — > e, 7. Nevertheless, as we will see, other processes 
like r — > /1, 7 and r — > e, 7 offer information complementary to that from /x — > e, 7, 
which in some cases can be extremely important. For r — > /1, 7, in particular, the 
present (future) upper limit is 1.1 x 10~ 6 (~ 10~ 9 ) [ p^JT9|1 . Finally, there exist interest- 
ing limits on other branching ratios, like BR(/i — > eee) < 10 -12 and BR(/i Ti —>■ e Ti) 
< 6 x 10 -13 |20"| J21f . The latter is also to improve in projected experiments p2 . These 



limits are certainly very small, but the theoretical predictions on the branching ratios 
are also much smaller due to the presence of extra electromagnetic vertices [p3[ 

From the results of the previous section, namely eq. (pOf) , we know that BR(/j — > 
oc \(Y+Y u )ij\ 2 . On the other hand, from the general equation (|HJ), it is clear that 
the matrix element (Y+Y„)jj will depend on several facts: the low-energy spectrum of 
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neutrinos (contained in D K ), the neutrino mixing angles (contained in U), the right- 
handed neutrino spectrum (contained in Dm) and the choice of the R matrix. Let us 
discuss them in order. 

The experimental (solar and atmospheric) data strongly suggest a hierarchy 
of neutrino mass-splittings, Ak 2 o( -C Aft 2 tm . Numerically, AK 2 tm ~ (1.4 — 6.1) x 
10~ 3 eV 2 /f| [|4[] (with central value around 3 x 10~ 3 eV 2 /t>|), while the value of Ak 2 o/ 
depends on the solution considered to explain the solar neutrino problem [p5|,p6j. The 



most favoured one from the recent analyses of data |24[ is the large-angle MSW solution 
(LAMSW), which requires A/« 2 oi ~ 3 x 10~ 5 eV 2 jv\. Other solutions are (in order of 
reliability from these analyses) the LOW MSW-solution, the small angle MSW-solution 
(SAMSW) and the vacuum oscillations solution (VO). They require, in eV 2 /f| units, 
Ak 2 oZ ~ 10~ 7 , 5 x 10~ 6 and 8 x 1CT 10 respectively. The corresponding uncertainties are 
discussed e.g. in ref. p3|. In any case, there are basically three types of neutrino spectra 



consistent with the hierarchy of mass-splittings ||27|| : hierarchical (nf <C /t 2 <C k^j, 



^ „, ^ ,^ 



"intermediate" (k\ ~ k 2 , ^> /t 2 ) and "degenerate" (k\ ~ k\ ~ k 2 ). We follow the 
usual notation in which k 3 corresponds to the most split mass eigenvalue, so that 

The hierarchical scenario seems perhaps the most natural one, since it resembles 
the other fermion spectra. Indeed, it is the natural spectrum in GUT theories, par- 
ticularly 50(10) models. In any case, it is the most extensively analyzed scenario 
in the literature and we will pay special attention to it. In this scenario one should 
take k 2 ~ Aft 2 tm , k 2 ~ Ak 2 ^, k\ ~ 0. The other two scenarios (intermediate and 
degenerate) are also interesting and we will address them. In particular, the degener- 
ate scenario is the only one compatible with a relevant cosmological role of neutrinos 
(provided n ~ C(l)eV/u|). 

Concerning the mixing angles, #23 and #13 are constrained by the atmospheric and 
CHOOZ data to be near maximal and minimal, respectively. The #12 angle depends 
on the solution considered to explain the solar neutrino problem. In the LAMSW case, 
612 is near maximal. In the other explanations of the solar neutrino problem #12 can 
be either minimal (SAMSW), or near maximal (VO and LOW). The precise preferred 
values for the angles depend on the sets of data used in the fits. For example, from 
the global analysis of ref.[|4[] these are tan 2 #13 ~ 0.005, tan 2 #23 ~ 1.4 and, for the 
LAMSW, tan 2 6*12 ~ 0.36. The uncertainties depend also on the type of fit, as well as 
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on the values used for the mass splittings. In summary, the two basic forms that U can 
present are either a single-maximal or (more plausibly) a bimaximal mixing matrix. 
Schematically, 
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(31) 



With regard to the spectrum of right-handed neutrinos, there are two basic possi- 
bilities: either they are hierarchical, i.e. M.\ <C M.% <C M.3, or they are degenerate, 
i.e. Ai\ ~ M. 2 ~ Both of them were considered in sect. 2, where it was also 

shown that the first possibility counts the case of an "inverse hierarchy" or any other 
ordering of the M. eigenvalues. 

Finally, R is the generic (complex) orthogonal matrix discussed in sect. 2. In 
principle, it contains three independent complex parameters, but we saw that this 
number can be usually reduced in a drastic way. 

Let us analyze now the predictions for BR(/j — > Ij,^) in the different scenarios 
above exposed. 

4.1 z/^'s and vrs completely hierarchical 
4.1.1 The generic case 

We consider first the case where both the left-handed and the right-handed neutrinos 
are completely hierarchical, i.e. 

D K ~ diag(0, « 2 , k 3 ), D M ^ diag(0, 0, M 3 ) . (32) 

We are neglecting for the moment contributions of 0(M.2). (Later on, we will discuss 
when those contributions would be relevant.) The general form of the Y v matrix in 
this case was already worked out in subsect. 2.1. If R32 7^ or R33 7^ (later on we 
will consider the case R32, R33 — 0), Y v is given by eq.fllED, so only the last row of Y v 
has sizeable entries. Thus, the relevant quantity (Y+Y„)^ has the form 



C^u)3iC^u)3j 



R llV^l U H 



1=2,3 



l'=2,3 



(33) 



Focusing on /1 — ► e, 7, and using the parametrization of R of eq.flTTD, eq.(^DJ) reads 
BR(/i^e, 7 )cx (Y+Y„) 2 r 

M3\C 2 \ 2 + C*V^^23] [siV^U* 2 + Ci^C/* 3 ] * (34) 



15 



Let us recall here that, since K2 and K3 represent the solar and atmospheric splitting, 
respectively, then K2 <C K3. Also, the matrix element C/13 ~ S13 <C 1, as #13 is the mixing 
angle constrained by CHOOZ. In the SAMSW scenario, U12 and U23 are also extremely 
small (see eq.@), so /x — > e, 7 is suppressed. However, in the other scenarios (LAMSW, 
VO and LOW) U is basically a bimaximal mixing matrix. Then U 12 ,U 23 ~ 1/V2, 
U22 — 1/2- Let us recall that this precisely includes the scenarios favoured from 
present data, i.e. LAMSW and LOW. 

The prefactor vVf 3 1 c 2 1 2 in eq . fl30"D can be written as 



M 3 \c 2 





Yo\ 


2 


1*1! 


2 K 2 + 


|ci 


^3 



(35) 



where |F | 2 is the largest eigenvalue of (Y+Y„). So, the right hand side of eq.([34|) 
depends just on two independent parameters, \Y \ and Q\ (the remaining parameters 
have a precise physical meaning and are measurable at low energy). If we further 
impose a unification condition, like |Y (Mx)| = |lt(Mx")|, BR(/x — > e, 7) depends just 
on one parameter, Q\. This conclusion could also be obtained from eq. (|18D . 

Now, as discussed in sect. 2, the results of ref. |10j correspond to take R — 1, i.e. 
Cj = 1. Then, BR(/i — >■ e, 7), as written in eq.(pl[) is suppressed by \Urs\ 2 <^ 1. The 
authors of ref. [|T(| considered however the contributions of 0(Ai 2 ) (which we have not 
considered yet). Then, working in the LAMSW scenario, they obtained a non- vanishing 
result for BR(/i — > e, 7), as a function of Ai 2 - What is apparent from eqs. (B3j, ^5] ) is 
that that result corresponds to a very special point in the parameter space, namely 
ci = 1. In general, ci 7^ 1, si 7^ 0, so the prediction for BR(/i — > e, 7) is much larger. 
For random values of 9\ we expect 

(Y+Y„) 21 ~ ^aiv^^23^ (36) 

This formula allows to compare the foreseeable predictions for BR(/i —>■ e, 7) in different 
scenarios, taking into account the different values for k%, Uy. Clearly, the largest 
values for BR(/i — > e, 7) are obtained for the LAMSW scenario (precisely, the most 
favoured one by the present analyses). For the LOW and SAMSW scenarios one expects 
predictions ~ 20 and ~ 1000 times smaller respectively. Indeed, for the LAMSW 
scenario the previous equation generically gives a branching ratio above the present 
experimental limits, at least for Yq = 0(1), as it occurs in the unified scenarios! 
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Figure 3: Branching ratio of the process /i — + e, 7 vs. the unknown angle 9\ for the case of hierarchical 
(left and right) neutrinos and two typical sets of supersymmetric parameters, as indicated in the 
plots. The dashed lines correspond to the present and forthcoming upper bounds. A top-neutrino 
"unification" condition has been used to fix the value of the largest neutrino Yukawa coupling at high 
energy. 9i is taken real for simplicity, so the two limits 9\ — 0, 7r of the horizontal axis represent the 
same physical point. In the upper plot the curves corresponding to to = 300 GeV and to = 500 
GeV are almost indistinguishable. The same is true in the lower plot for the curves corresponding to 
mo = 100 GeV, to = 300 GeV and m = 500 GeV. 
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This is apparent from the plots of Fig. 3, where it is shown BR(/i — > e, 7) vs. 9\ 
in the LAMSW scenario for different typical sets of supersymmetric parameters, using 
a unification condition for the largest neutrino Yukawa coupling, i.e. |lo(^x)| 2 = 
\Y t (Mx)\ 2 , with Mx = 2 x 10 16 GeV. Incidentally, the value of the supersymmetric \i- 
parameter is extracted (here and throughout the paper) from the electroweak breaking 
condition, choosing /i > 0. We have not shown the results for < since they are 
always very similar. From Fig. 3, it is clear that most of the parameter space (i.e. the 
range of #1) is already excluded by the present bounds on BR(/i — > e,j), or it will 
become probed by the next generation of fi — > e, 7 experiments (PSI), scheduled for 
2003. For the sake of clarity, we took Q\ real in Fig. 3. The results for 9\ complex 
are very similar, as can be seen e.g. from eq. (|36|) . Fig. 4 is completely analogous to 
Fig. 3, but for tan/? = 10 (instead of tan/3 = 3). It illustrates the strong dependence of 
BR(/i — > e, 7) on tan/3 in the expected way (~ tan 2 /3). An alternative representation 
of the allowed parameter space is given in Fig.5, where Q\ is given in logarithmic 
units (the set of supersymmetric parameters corresponds to that of the second plot of 
Fig. 3). These can be more meaningful in particular model-building constructions. Of 
course, the basic results remain: 6\ should be < 0.12 (< 3 x 10~ 3 ) with the present 
(future) upper bounds on BR(/i — > e, 7). This corresponds to Y31 < 0.04Y 33 (Y31 < 
9 x 10~ 4 Y 33 ). This gives a measure of the constraints on the textures. (In a logarithmic 
plot it is not possible to say what percent of the parameter space is allowed since this 
includes the origin.) 

It is worth- noticing from the curves of Figs. 3, 4, the existence of two values of §1, 
namely 9\ = 0, 9\ ~ 2, in whose neighborhood the value of BR(/i — > e, 7) is drastically 
suppressed. We will comment shortly on this feature. For large enough values of 
the universal soft scalar mass, m , BR(/z — > e, 7) decreases, as expected: for large 
m , the branching ratio goes as ~ "^o 4 , as it is apparent from Fig. 6. Typically, for 
tan /3 = 3 BR(/i — > e, 7) does not fall below the present or forthcoming experimental 
upper bounds until m ^ 1.5 TeV, i.e. beyond the reasonable range to avoid fine- 
tuning problems. For larger tan (3, this extremal value increases, since BR(/i — > e, 7) ~ 
tan 2 f3/rriQ. Notice that there is a very narrow window at around 120 GeV, for A = 
GeV, where the branching ratio is very suppressed. This happens for every value of §1 
and is the result of a cancellation between the chargino and the neutralino amplitudes. 
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Figure 5: The same as Fig. 3 (lower plot) using a logarithmic scale for d\. 

Finally, if one relaxes the unification condition on |YqCMx:)I; BR(/i — > e, 7) logically 
decreases, as shown in Fig. 7. Anyway, it is noticeable that even for Y^iMx) ten times 
smaller than Y t (Mx), the conclusion remains that most of the parameter space will be 
probed in the forthcoming experiments. 

4.1.2 Special textures 

Let us now turn to the two values of 6\ where BR(/i — > e, 7) is suppressed. They 
correspond to the values for which (Y+Y i/ ) 2 i, as given by eq. (|3"3"D , is small. Notice 
that there are two possible ways to fulfill this requirement, namely (Y„)3i ~ or 
(Y^)32 — 0, which translate into 



tan#i 
tan#! 



K 2 U{ 2 



K 2 V{ 2 
K 2 V 2 \ 







(37) 
(38) 



K 2 [/| 2 

where we have used eqs.(|18|, |8|). We do not know any reason why any of the previous 
conditions should be satisfied in a particular model, although of course it is not ex- 
cluded. To this respect, the first condition, i.e. Q\ ~ 0, seems to be more "natural" 
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Figure 6: Branching ratio of the process \i — > e, 7 vs. the universal scalar mass, mo, in the case of 
hierarchical (left and right) neutrinos, and for the two sets of supersymmetric parameters of Fig. 3 
(with tan/3 = 3) and 9\ = 0.5. The dashed lines correspond to the present and forthcoming upper 
bounds. Except for a narrow range at low energy, the prediction does not fall below the present upper 
bound until m > 1.5 TeV. The same occurs for other values of 9\, except for the two special values 
seeable in Fig. 3 and discussed in subsect. 4.1.2. 
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A =250 GeV, M =250 GeV, m =100 GeV, tan(3=3, ^>0 



10" 



10" 



?*- 

CD 

A 

l_10" 



10" 



1 I 1 I 1 I 

Y ~Y 

\ x ^ t 


l Y-Y/5 \ 








lit, 





10" 

12 3 

«, 

Figure 7: The same as Fig. 3 (lower plot) for different values of the largest neutrino Yukawa coupling, 
Yq, at the "unification" scale, Mx- Y t denotes the value of the top Yukawa coupling at that scale. 



to obtain than the second one: note that 9i = 0, i.e. R = 1, corresponds to the case 
where all the mixing can be attributed to the sector of charged leptons, as discussed in 
sect. 2, which is a physical possibility. It is interesting to explicitly show the textures 
of Y u that correspond to these 'privileged' values of 9\ from eqs. (|3T|, |38|). They read, 
respectively 



(39) 
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Here we have used Iml/ii 



numerical values of the entries in (|39 



(40) 



and S13 1, so that ImV2i,ImV3i 1. Besides, the 
correspond to a bimaximal mixing V matrix, 



which includes the preferred LAMSW scenario. Let us stress that the previous textures 
are to be understood in the basis we have chosen to work, i.e. where both Y e and A4 
are diagonal, and the latter has the eigenvalues ordered in the usual hierarchical way 
shown in eq. (|32|) . It is interesting to note that texture (|39| ) has been advocated in 
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ref.[0 on different grounds. 

Some comments are in order here. As was discussed in subsect. 2.1, the first two 
rows of Y„, though suppressed, are normally important to reconstruct k from Y„, M.. 
In particular, for texture fl40|) the precise values of the first two rows are crucial to 
reproduce the bimaximal structure of k. On the other hand, this is not necessarily the 
case for texture ([39]). Actually, it is remarkable that texture fl39|) works fine, both to 
reproduce k and to give suppressed BR(/i — > e, 7), even if it is given in a basis where 
M. is not diagonal. 

In general, it is clear from eq.(^) that the k — > process is suppressed at the 
points in the parameter space where (Yj,^ ~ or (Y v ) 3 j — 0. Therefore, there are 
three special textures: the two given in the above eqs.(|3^, [5]), which correspond to 
(Yy)3i ~ and (Y y ) 32 ~ 0, and the one corresponding to (Y„) 33 ~ 0. The latter 
implies 

tan(il H = /E£ , (4i) 

and the associated texture is 

/ 0\ / 0\ 
Y^oc oc 00 , (42) 
V^2i -V n 0j [-1 y/2 Oj 

where, again, the numerical values values correspond to the bimaximal scenario. If the 

texture of Y u is one of the special forms (|39| , ^Oj ^2]) , there are two U — > Zj, 7 processes 

suppressed and one unsuppressed. This means that for those privileged cases where 

li — » e, 7 is unusually suppressed, the r — > 7 and r — > e, 7 processes are relevant 

and may give the most important restrictions on models. In particular, if Y u is as in 

eq.© or eq.flP (i.e. (Y„) 31 ~ or (Y,) 32 ~ 0), then BR(r -> /i, 7 ) or BR(r -> e, 7 ) 

respectively are unsuppressed. 

We have illustrated this interesting aspect by plotting BR(r — > 7) vs. ^1 in Fig. 8 

for one of the typical supersymmetric scenarios considered in Fig. 3. Notice that for 

the special value Q\ ~ 0, at which BR(/i — > e, 7) is suppressed, BR(r — > /i, 7) is not, 

lying above the forthcoming experimental upper bound. 

4.1.3 The contribution oc A4 2 

So far, all the results of this subsection have been obtained neglecting the value of M.2 
(see eq.(B2|)). Let us consider here a value of M. 2 7^ 0, though smaller than A4 3 , as it 
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Figure 8: The same as Fig. 3 (lower plot), but for the branching ratio of the process r — ► /x, 7. The 
dashed lines correspond to the present and forthcoming upper bounds. 



was done in ref . [[Kj . 

Clearly, for generic values 9 i: the contributions already computed [see eqs.(]3"4]- 
P^)] will be the dominant ones. However, if 9\ corresponds to one of the particular 
values of eqs.(|3"7|, |3"B"D , the contributions ~ M. 2 can be relevant or even dominant. 
Normally, this will happen only for §1 quite close to eqs.fl37j, |38|) . Suppose for example 
that R ~ 1, i.e. s, small and real. This corresponds to eq. (|3"TD and it is probably 
the most interesting case (and the one considered in ref. ||10|| ). Then (Y+Y„)y = 
(Y l/ )* 2i (Y u ) 2 j + (Y i/ )3 i (Y iy ) 3i with Y u given by eq.([L3|). Keeping only terms up to 0(s») 



we get 



(Y+Y^i ~ M 2 k 2 U* 12 U 22 + MzKjjyJn 

+ h{M s - M2)y/w*{TjyJn + U* 2 U 23 ) , 



(43) 



so that in a bimaximal case, even taking U13 = 0, the part proportional to will be 



indeed the dominant one except for s% < (A^/A^Xy ^2/^3), i-e. very close to zero. 
For the specific case i? = 1, i.e. Si = 0, we get (Y+Y„) 2 i — A4 2 k 2 U^ 2 U 22 . Taking 



into account from eq.( 



|Yo| 2 = the value of (Y+Y !/ ) 2 i is exactly as the one 

24 



we will obtain in another context in subsect. 4.2, see eq.(^) below, but multiplied by 
.M2/.M3. In consequence, the plots of BR(/i — > e, 7) are identical to those of Fig. 9 
below, but multiplied by (.M2/.M3) 2 . For M. 2 j M.^ > 10~ 2 the branching ratio is still 
testable in the forthcoming generation of experiments. 

4.1.4 i? 32 ^ R33 ^ 

So far we have assumed -R32 7^ or R 33 7^ in eq. flT3|) . Recall that R is an unknown 
orthogonal matrix. The experimental data do not restrict its form, although R may 
be much more definite in specific models. We consider now the special case when 
i? 32 ~ -R33 ~ 0. It should be kept in mind that this represents a narrow region of the 
whole parameter space of R. 

However, as we are about to see, this possibility includes the particular case where 
all the leptonic flavour violation can be attributed to the sector of right-handed neu- 
trinos. This happens when there exists a basis in which Y e and Y u get simultaneously 
diagonal, while A4 gets non-diagonal. This is physically interesting and it can occur 
in particular models. As explained in sect. 2, in the basis we have chosen to work (Y e 
and A4 diagonal), this implies Y v = WDy, where Dy is a diagonal matrix and W is a 
unitary matrix. Hence, Y u + Y u = D\ Y 2\- Then, from the general equation (|TB| ) or (|H|), 
we get 

R+D M R = D^-jU+D^UD^-t . (44) 

If U is of the bimaximal type, as it occurs in the preferred LAMSW scenario, then 
Un 7^ for any i [see eq. (pi])]. This means that, in first approximation, the only non- 
vanishing entry of the right hand side of eq.([P]) is the (1,1), so the same holds for the left 
hand side. Since Dm is hierarchical, (R + D^R)^ — R%iM-zRzj ■ Thus, -R32 — -R33 — 
q.e.d. A special case occurs when D\y*\ = |^q 2 |1- Then eq.(f£|) becomes R + D M R = 
|y o 2 |-D K -i for any form of U; and R turns out to be simply the "permutation matrix" 
that interchanges the 1- and 3-planes: R = R p = {(0,0, ±1), (0, ±1,0), (±1,0,0)}. 
Of course, in all these cases, since Y^ + Yj, is diagonal, there is no RG generation of 
non-diagonal soft terms, and the predictions for BR(Zj — > ^,7) are as in the SM, i.e. 
negligible. 

Nevertheless, the possibility R 32 — R33 — includes other cases different than 
the previous (almost trivial) one. As mentioned, this represents a small region of the 
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parameter space, which in principle has no special physical significance; but for the 
sake of completeness we will analyze it. R will admit the parametrization 



R ~ ts c , (45) 




where 9 is an arbitrary complex angle and the sign of the ±1 entry is not correlated to 
the signs appearing in the second column. Thus, from eq. (|T3|) , and keeping only the 
dominant terms, 



±y / Mic y fR^ VMiSy/t^' 
Y v ~ I T^FM~ 2 sJTT 2 ^[M~2cJ^i I U + , (46) 



±y/M^i 

Therefore, the values of (Y u + Y,y)ij, i ^ j, which are the relevant quantities for U — > 
lj, 7, depend on the value of s. In general, if a unification condition is imposed for the 
largest eigenvalue of Y u + Y u , then (Y u + Y u ) 2 i is quite large and BR(/i — > e, 7) is above 
the experimental limit. 

Let us illustrate the last point with an example: suppose that s = ±1 in eq. fl4*6]) . 
Then the eigenvalues of Y u + Y u are {M.3K1, M-i^}- If they are degenerate, 

Y u + Y u is diagonal and we are in the above-considered case: all the mixing can 
be attributed to the sector of right-handed neutrinos. If not, since U13 ~ 0, then 
(Y u + Y u ) 2 i — + C^i2^22-A42«2- This is at least as large as the matrix 

element computed in the generic case, eq . (|36|) , thus leading to unacceptably large 
BR(/i — > e, 7). For generic values of s the analysis is more involved, but at the end of 
the day the conclusion is similar. 

4.2 z/x's hierarchical and vrs degenerate 

As was discussed in subsect. 2.2, if the R matrix is real, this scenario is extremely 
predictive. In particular Y+Y„ is given by eq.(f2~l~l), which does not depend on the form 
of R. More precisely 

{Y+YJy = Mk 1 U u U+ ~ M [k 2 U i2 U* 2 + k 3 U i3 U* 3 ] , (47) 

where k 2 = k so i and k 3 = K atm . The free parameter M. can be absorbed in the value 
of the largest eigenvalue of the Y+Y„ matrix, |y | 2 — M.K3- In particular, the matrix 
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element (Y+Y I/ ) 2 i becomes 

(Y+Y„) 2 i ~ Mk 2 U 22 U* 2 = \Y \ 2 —U 22 U* 2 , (48) 
where we have used U13 ~ 0. Comparing (^) with (|36l), we see that in this scenario 



we expect predictions for BR(/i — > e, 7) and the other processes of the same order as 
for hierarchical neutrinos (more precisely, ~ K 2 j k 3 times smaller), but with no special 
textures where the branching ratio becomes suppressed. All this is illustrated in Fig. 9. 
for the LAMSW scenario using |Y (M X )| = \Y t {M x )\. 

From the plots, it is clear that in this scenario BR(/i — > e, 7) is already above the 
present experimental limits except for a rather small region of m -values which should 
be probed by the next generation of experiments. 

If R is complex, the analysis is more involved since it contains more arbitrary 
parameters; but in general the conclusion is the same: BR(/z — > e, 7) is at least of the 
same order as in the real case. To see this note that from the general equation (|i~4|) 



Y+YJ01 = M 



1<2 



M U 23 U* 12V ^R* q3 R q2 + U 22 U* 12 K 2 R* q2 R q2 . (49) 



Normally, this is at least of the same order as in the real case, eq. fl4"8"|) , as we have 
checked numerically. Now there exists, however, the possibility of a (fine-tuned) can- 
cellation between the various contributions of (£|^). In particular such cancellation will 
occur if all the mixing can be attributed to the sector of right-handed neutrinos, which 
implies Y+Y u diagonal. 

4.3 vls quasi- degenerate 

As it was discussed in subsect. 2.3, in this case it is logical to assume that M. has 



degenerate eigenvalues. If, furthermore, the R matrix is real, we get from eq. (|T4]) 



Y+Y, = MUD K U + , (50) 

Incidentally, for R real the requirement that for quasi-degenerate neutrinos the two CP 
phases, and 0', are opposite in order to be consistent with i/-less double (3 decay p8| , 
implies that U and thus Y u is complex, see eqs.(|8], |9|, |13| ). 

Since Ki ~ k 2 ~ k 3 = k (typically k = C(leV /v%)), the eigenvalues of Y+Y^ 
are quasi-degenerate ~ Ain = \Y \ 2 . Clearly, eg-floPP presents a GIM-like suppression, 
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M =250 GeV, tan(3=3, n>0 
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10" 
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m 

Figure 9: Branching ratio of the process /i — ► e, 7 vs. the universal scalar mass, too, for the case 
of hierarchical (degenerate) left (right) neutrinos, R real (see subsect. 4.2) and two typical sets of 
supersymmctric parameters, as indicated in the plots. The dashed lines correspond to the present and 
forthcoming upper bounds. A top-neutrino "unification" condition has been used to fix the value of 
the largest neutrino Yukawa coupling at high energy. The two plots illustrate two different values of 
tan/3. The curves do not fall below the present bound until m > 1.6 TeV (uper plot) and m > 2.9 
TeV (lower plot). 
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so the off-diagonal entries of Y+Y„ are small and proportional to the differences of 
neutrinos masses. More precisely, for i ^ j 

(Y+Y„)y = MU ip U* p k p = MU ip U* p (K p - « 3 ) . (51) 

In particular, the matrix element (Y+Y jy ) 21 , which is the relevant one for \x — > e, 7, 
becomes 

(Y+Y u ) 21 =M[U 21 U* u (k 1 -k 3 ) + U 22 U* 2 (k 2 -k 3 )} . (52) 
Since U\ 3 — 0, [] from unitarity follows that U^U^ ~ —U 22 U{ 2 . So 

(Y+Y„) 21 = M [U 21 U* 1x {k x - k 2 )\ = IKoft^tfn^k . (53) 



K 2 



The factor A ^f°' is at most ~ 10~ 5 (this occurs for the LAMSW solution to the so- 
lar neutrino problem), which implies a drastic suppression of (Y+Y„) 2 i, and thus of 
BR(/i — > e,7) in this kind of scenario. More quantitatively, eq.(|53"D is identical to 
eq.(f48|) (i.e. the scenario where Ul's are hierarchical and vr's degenerate), multiplied 



by k~ 2 y AK 2 ol AK\ tm . This is a factor ~ 10~ 4 for the LAMSW, which means that all 
the plots representing BR(/i — > e, 7) in that scenario (Fig. 9) are identical here, but 
with the vertical axis re-scaled eight orders of magnitude smaller. 

If the R matrix is complex, then R + R is in general non-diagonal, so the Y+Yj, 
matrix becomes 

Y+Y„ = MUD^R + RD^U+ , (54) 

which may have sizeable off-diagonal entries. Hence, even for quasi-degenerate neutri- 
nos, (Y+Y v ) 2 i, and thus BR(/i — > e, 7), could be very large if R is complex. 

Finally, we would like to comment on the scenario where the neutrinos are partially 
degenerate, i.e. when their masses present an inverse hierarchy: k 3 <C K\ ~ k 2 = k. [It 
is convenient to maintain k 3 as the most split mass eigenvalue, so that the interpretation 
of #23, #12 a s the atmospheric and solar angles respectively remains valid.] Notice that 
in this case k 2 ~ AK 2 tm . 



^Since we are working here with quasi-degenerate neutrinos, one may wonder whether the low- 



energy condition U13 ~ might change at the M. scale. However, at is was shown in rcf. [29 , the 
U13 element is stable under the RG running between M. and low-energy for any phcnomcnologically 
viable scenario. Therefore, the argument that follows is completely general. 
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Now, the discussion of the size of (Y+Y Jy ) 2 i, and thus of BR(/i — > e, 7), is similar 
to that in the complete degenerate case, though the conclusions are different. In 
particular, suppose that M. has degenerate eigenvalues, Dm = diag(.M, A4, M.), and 
that the R matrix is real. Then eq. (|50|) holds, so Y^Y U has two large quasi-degenerate 
eigenvalues ~ M.K = \Y \ 2 . Likewise, all the discussion between the eqs. (|5l"r) -(|53"|) 
and the equations themselves remain valid. Therefore, the matrix element (Y+Y^i 
becomes 

(Y+Y„) 21 = \Y \ 2 U 21 U* n ^ = \Y \ 2 U 21 U^^^- . (55) 

K ^ K atm 

Again, this equation is identical to eq.fl48|), multiplied now by 4 / f K £ oL . This suppression 

V atm 

factor is not as strong as before. Namely, it is ~ lO^ 1 for the LAMSW, which means 
that all the plots representing BR(/i — ► e, 7) in this scenario are as in the one where 
ul's are hierarchical and ur's degenerate, re-scaled by a factor ~ 10 -2 . Hence, from 
Fig. 9 it is clear that BR(/i — > e, 7) for the partially degenerate case should be testable 
within the next generation of experiments. 

The matrix element (Y+Y^) 31 gets a similar suppression. However (Y+Yj,^ is 
not suppressed at all: 

(Y+Y,) 32 = M [UsiU^Kt - k 3 ) + U 32 U; 2 (k 2 - k 3 )} ~ Mk [U 3l U* 21 + U 32 U* 2 ] ,(56) 

where M.K = \Yq\ 2 . This is a similar size to that from the hierarchical scenarios. So, 
for the partially degenerate case, BR(r — > fi, 7) should also be testable in the next 
generation of experiments. 

Again, if the R matrix is complex Y+Y^ is given by eq. (|54D and we expect much 
larger values for (Y+Y i/ ) 21 , (Y+Y l/ ) 31 . 

5 Summary of predictions for BR(/^ — > lj, 7) 

If the origin of the neutrino masses is a see-saw mechanism, implemented in a super- 
symmetric theory, the soft breaking terms get off-diagonal contributions through the 
RG running, even if one starts with a universality condition for the soft terms. E.g. 
for the slepton doublet mass-matrix, at the leading-log approximation, 

H). * 8^ (3m ° + A o)(Y, + Y,),- log ^ , (57) 
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where Y u is the matrix of neutrino Yukawa couplings, M is the typical (Majorana) 
mass of the right-handed neutrinos and Mx is the initial scale at which universality is 
imposed. Therefore, the predictions on BR(/j — > ^,7) are directly linked to the size 
of (Y+Y„)ij. Moreover, the larger (smaller) the initial scale at which universality is 
imposed the larger the branching ratios. Also, due to the structure of the diagrams, 
BR(/j — > has a strong dependence on tan/3 (~ tan 2 /?). The most general form 
of Y„ and thus of Y+Y u , in terms of measurable low-energy parameters has been 
determined in sect. 2. In particular, working in the flavour basis in which the gauge 
interactions and the charged-lepton Yukawa matrix, Y e , are flavour-diagonal, 

Y+Y, = UD^R+D M RD^U + (58) 

Here is the diagonal matrix of positive neutrino mass eigenvalues (k = M. u / (H®) 2 ), 
while U contains the mixing angles and CP phases [see eqs.(|| [7p]. In addition, 
Y+Y^ depends on unknown parameters: the three positive mass eigenvalues of the 
righthanded neutrinos (contained in Dm) and the three complex parameters defining 
the complex orthogonal matrix R. In practical cases, however, the number of relevant 
free parameters becomes drastically reduced. 

Next, we enumerate the possible low-energy scenarios and the corresponding pre- 
dictions for (Y+Yj,)y and BR(I, — > lj, 7), as has been worked out in detail in sects. 2, 
3, 4. 

z/^'s and z/r's completely hierarchical 

• If R is a generic matrix, (Y+Y^)^ is given by eq. (|3"B|) . Using the parametrization 
of R given in eq . (|i~7D , which is sufficiently general for this case, one obtains in 
particular 



(Y+Y,) 21 ~ 1 cUhVwGUnUi* (59) 

|Sl| «2 + \Cl\ Z K 3 

Here |Y | 2 is the largest eigenvalue of Y+Y„ [see eq. (|35|)] and 6\ is an arbitrary 



complex angle. 

For the LAMSW scenario the previous equation generically gives BR(/i — > e, 7) 
above the present experimental limits, at least for Y Q = 0(1), as it occurs in the 
unified scenarios. This is illustrated in Figs. 3, 4, 5, 6 using |1q(Mx)| = |F((Mx)|, 



31 



with M x = 2 x 10 16 GeV. But even for Y = (9(10 _1 ), most of the parameter 
space will be probed in the forthcoming generation of experiments experiments 
(see Fig. 7). 

• The only exceptions to the previous result are 

— If 61 takes one of the two values given in eqs. (|37| , |3~8l), which correspond 
to the textures of eqs. ([39], f40|) . Then, BR(/i — > e, 7) is drastically reduced, 
but other processes, as BR(r — > /i,7), are not, normally lying above the 
forthcoming experimental upper bound, see Fig. 8. 

— If R is such that Y+Y^ is diagonal. This means that in a certain basis 
Y e and Y u are diagonal whereas the right-handed Majorana mass matrix, 
Ai, is not. This requires a very special form of R, which in particular has 
R32, R33 — 0. Of course, in this special case there is no RG generation of 
non-diagonal soft terms, and the predictions for BR(/j — > /,,-, 7) are as in the 
SM, i.e. negligible. 

z/^'s hierarchical and z/r's degenerate 

• If R is real, this scenario is very predictive. (Y+Yj,)^- is given by eq.(^). In 
particular, taking into account U13 ~ 0, 

(Y+Y„) 21 ~ Mk 2 U 22 U* 2 = \Y \ 2 —U 22 U* 2 , (60) 

^3 

where |Y" | 2 — M-^3 is the largest eigenvalue of Y+Y^. 

The corresponding BR(/i — > e, 7) is shown in Fig. 9. for the LAMSW scenario, 
and using = |¥t(M_x-)|. [] The branching ratio turns out to be already 

above the present experimental limits except for a rather small region of m 
values which should be probed by the next generation of experiments. Here are 
not special textures where the branching ratio becomes suppressed. 

• If R is complex, the analysis is more involved since it contains more arbitrary 



parameters; (Y+Y^) 2 i is given by eq. fl49|) . But in general the conclusion is the 



same: BR(/x — > e, 7) is at least of the same order as in the real case. Now there 



7 For other values of |Yo| note that BR(^i — > e, 7)) ~ |^o| 4 - Actually, this dependence is milder, 
since the larger |5^>(Mx-)|, the faster it decreases with the scale. 
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exists, however, the possibility of a (fine-tuned) cancellation between the various 
contributions of (fj5|). In particular such cancellation will occur if all the mixing 
can be attributed to the sector of right-handed neutrinos, which implies Y+Yj, 
diagonal. 

vls quasi- degenerate 

• If R is real, (Y+Y„)ij is given by eq. (|5H|) . In particular (Y+Y v ) 2 i is given by 

(Y+Y„) 21 = M [U 21 UtM - k 2 )] = \Y \ 2 U 21 U^^- . (61) 
where, again, | Yq | 2 ~ A4k is the largest eigenvalue of Y+Yj,. This equation is 



identical to eq. (|60D , multiplied by k~ 2 J Ak^Ak^. This is a factor ~ 10~ 4 for 
the LAMSW. Therefore all the plots representing BR(/i — > e, 7) in the previous 
scenario (Fig. 9) are valid here, but with the vertical axis re-scaled eight orders 
of magnitude smaller. Consequently, BR(/i — > e, 7) is naturally suppressed below 
the present (and even forthcoming) limits. 



If R is complex, Y+Y^ is given by eq. fl54|) , which may have sizeable off-diagonal 



entries. Hence, BR(/i — > e, 7), could be very large in this case. 
If the (quasi-) degeneracy is only partial: k 3 < K! ~ k 2 e k ~ \I An, 



.2 

(Y+Y I/ ) 2 i is given (for R real) by eq. (|55|) . Again, this equation is identical to 
eq. (pCf) , multiplied now by J f V' ■ This represents a suppression factor ~ 1CT 1 

y atm 

for the LAMSW, which means that all the plots representing BR(/i — > e, 7) are as 
those of Fig. 9, re-scaled by a factor ~ 10~ 2 . As a consequence, BR(/i — > e, 7) for 
this partially degenerate scenario should be testable within the next generation 
of experiments. The conclusion is similar for BR(r — > /i,7). 

For generic complex R, the value of BR(/i — > e, 7) does not get any suppression 
and falls naturally above the present experimental limits. 

6 Concluding remarks 

If the origin of the neutrino masses is a supersymmetric see-saw, which is probably 
the most attractive scenario to explain their smallness, then the leptonic soft breaking 
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terms acquire off-diagonal contributions through the RG running, which drive non- 
vanishing BR(/j — > lj, 7). These contributions are proportional to (Y+Yj,)^-, where Y„ 
is the neutrino Yukawa matrix, 

Therefore, in order to make predictions for these branching ratios, one has first 
to determine the most general form of Y„ and Y+Y„, compatible with all the phe- 
nomenological requirements. This has been the first task of this paper, and the result 



is summarized in eqs.(|13|, |14|) 



Then, we have shown that the predictions for BR(/i — > e, 7) are normally above the 
present experimental limits if the three following conditions occur 

1. The solution to the solar neutrino problem is the LAMSW, as favoured by the 
most recent analyses. 

2. Y (M X ) = where |F | 2 is the largest eigenvalue of (Y+Y^). This occurs 
e.g. in most grand-unified scenarios. 

3. The soft-breaking terms are generated at a high-energy scale, e.g. Mx, above 
the Majorana mass of the right-handed neutrinos, M. 

These conditions are very plausible. In our opinion, the most natural scenarios fulfill 
them, but certainly there exists other possibilities. E.g. concerning condition 1, the 
solution to the solar neutrino problem could be LOW. In this case, BR(/i — > e, 7) is 
(~ 20 times) smaller, but still testable in the forthcoming experiments. Condition 
2 is satisfied in SO(10) models, but it could be unfulfilled in SU(5) models or in 
string scenarios with no GUT group. Still, we have noted in this work that even 
for Y = (^(lCr 1 ), most of the parameter space will be probed in the forthcoming 
generation of experiments (year 2003). Finally, condition 3 seems a natural one, as 
the see-saw mechanism requires the existence of high-energy scales, and it is hardly 
compatible e.g. with models where the fundamental scale is O (TeV). Actually, if the 
initial scale at which universality is imposed is larger than the ordinary GUT Mx scale 
(as it is common in gravity mediated supersymmetry breaking), the predictions for the 
branching ratios increase. However, it may happen that supersymmetry is broken at a 
scale below M. This is the case of gauge-mediated scenarios, where there would be no 
generation of off-diagonal leptonic soft terms through the RG running. In this sense, 
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the gauge-mediated models are free from /i — ► e, 7 constraints and represent a way-out 
to the previous conclusion. 

Even under the previous 1-3 conditions, there are physical scenarios compatible with 
the present BR(/x — > e, 7) experimental limits. Namely 

• Whenever all the leptonic flavour violation can be attributed to the sector of 
right-handed neutrinos. This happens when there exists a basis in which the 
gauge couplings, Y e and Y u get simultaneously flavour-diagonal, while the right- 
handed Majorana matrix, A4, gets non-diagonal. In this case there is no RG 
generation of non-diagonal soft terms, and the predictions for BR(7j — ► ^,7) are 
as in the SM, i.e. negligible. 

This possibility certainly exists and might be considered an attractive one. How- 
ever, it should be noticed that it does not occur in the quark sector since the Y u 
and Yd Yukawa matrices cannot be simultaneously diagonalized in a basis where 
the gauge couplings remain diagonal. In our opinion, it would be a big chance 
that this happens in the leptonic sector. 

• In the scenario where both, the left-handed and the right-handed neutrinos have 
hierarchical masses, there are two special textures, shown in eqs. fl3"P|, fTO|), where 
BR(/i — > e, 7) becomes suppressed. The first one is probably the most appealing, 
since it takes place for the simple case R — 1 in eq.([H|). This is equivalent to 
assume that there exists a basis of L, and Vm in which ~Y V and Ai are simultane- 
ously diagonal, though not Y e . In other words, all the leptonic flavour violation 
can be attributed to the sector of charged leptons. 

Then, we have shown that for this privileged texture BR(r — > fi, 7) is not sup- 
pressed, lying in general above the forthcoming experimental upper bound. 

• If the left-handed neutrinos are quasi-degenerate, then BR(Z; — > are natu- 
rally suppressed. Namely, if the R matrix in eq.(|T3|) is real, there is a GIM-like 
mechanism that drives BR(7j — > below the experimental limits, even the 
forthcoming ones. (This is not the case, however, if R is complex.) 

If the degeneracy is partial (case of inverse hierarchy of neutrino masses), BR(/i — ► 
e, 7) is somewhat suppressed, but still testable in the forthcoming experiments. 
The latter is also true for BR(r — > fi, 7) 
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In our opinion, if the above 1-3 conditions are fulfilled, the scenario of quasi- 
degenerate neutrinos and the one with gauge mediated supersymmetry breaking repre- 
sent the most plausible explanations to the absence of fi — > e, 7 observations, specially 
if the absence persists after the next generation of experiments. 

As a final conclusion, the discovery of neutrino oscillations makes much more plausi- 
ble the possibility of observing lepton-flavour- violation processes, specially // — > e, 7, if 
the theory is supersymmetric and the neutrino masses are generated by a see-saw mech- 
anism. Large regions of the parameter space are already excluded on these grounds, 
and there exists great chances to observe fi — ► e, 7 in the near future (PSI, 2003). This 
means that, hopefully, we will have signals of supersymmetry before LHC. 
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Appendix 



In the appendix, we summarize the relevant formulas concerning RGEs and the coef- 
ficients of the amplitude appearing in eq. fl2"T|) . 

Between Mx and M the evolution of the Yukawa couplings is given by 



dY v 
dt 



16tt 2 



3gi + -g( - T 2 Is - 3Y + Y„ + Y e + Y 



(62) 



dY e 



3g 2 2 + 9 -g\ - T x ) I 3 



Y+y 



3Y+Y e 



(63) 



dY 



u 



dt 



16vr 2 



-Y 



u 



13 
15 



9 2 i + Zgl + y 2 3 2 - 2a) I 3 - 3Y+Yu - Y+ Y D 



(64) 



where 



T! = Tr(3Y+Y D + Y+Y e ), T 2 = Tr(3Y+Yu + Y+Y, 



(65) 
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and 



dM 



M(Y V Y U + ) T + Y U Y U + M 



dt 8tt 2 

The renormalization group equations for the soft masses and trilinear terms are: 
dmj _L_ 
dt 



(66) 



dm 2 

dt 



dm 2 

dt 
dm 2 H2 



dt 



dA e 



dA v 



dA 



u 



dt 



where 



16tt 2 



(ml Y+Y e + Y+Y e m 2 ) + (m 2 Y+Y, + Y+Y.m 2 ) 



+2 (Y+m2Y e + Y+Y e + A+A e ) + 2 (Y+mX + m^Y+Y„ + A+A„) 



^2|ji / j-|2 iC 2|ji / j-|2\ t 3 2 , 



- ( -g{ \Mtf + 6^ |M 2 | Z J I 3 - -^Sl 3 



16tt 2 
24 

~5 
1 



2 (m 2 Y e Y+ + Y e Y+m e 2 ) + 4 (Y e m 2 Y+ + m 2 Hl Y e Y+ + A e A+) 



^ 1 2 |M 1 | 2 I 3 + ^I 3 



(67) 



(68) 



16tt 2 
1 

16tt 2 



2 (m 2 Y,Y+ + Y,Y+m 2 ) + 4 (Y,m 2 Y+ + m 2 H2 Y v Y+ + A„A+)] ,(69) 
6Tr (Y+(mJ + + m^I 3 )Y u + A+ A v ) 



+2Tr (Y+(m 2 + m 2 + m| 2 I 3 )Y„ + A+A„) 
-(^ 2 |M 1 | 2 + 6y 2 |M 2 | 2 )+^", 

^ [{-^i 2 - 3^ 2 2 + 3Tr(Y+Y D ) + Tr(Y+Y e )} A e 

+2 {- jj^Mi - 3^ 2 2 M 2 + 3Tr(Y+ A D ) + Tr(Y+A e )} Y e 
+4(Y e Y+A e ) + 5(A e Y+Y e ) + 2(Y e Y+A,) + (A e Y+Y,) 
[|-^ 2 - 3g 2 2 + 3Tr(Y+Yu) + Tr(Y+Y,)} A u 

+2 j-^Mi - 3<? 2 M 2 + 3Tr(Y+Au) + Tr(Y+A,)| Y, 

+4(Y,Y+A„) + 5(A„Y+Y V ) + 2(Y„Y e + A e ) + (A„Y e + Y e )" , 

16 1 
-9l ~ Zgl - -gl + 3Tr(Y+Yu) + Tr(Y+Y„) j Au 

16 



(70) 



(71) 



(72) 



16tt 2 



13 
15* 



+2 {--^Mi - 3s 2 M 2 - -g 2 3 M 3 + 3Tr(Y+A u ) + Tr(Y+A,)} Y v 
+4(Y U Y+A U ) + 5(A U Y+Y U ) + 2(Y Y+ A D ) + (AuY+Y D ) 



S = Tr(mJ + - 2m 2 - ml + m 2 ) - m 2 H + mj 



(73) 



(74) 



Here # 2 and g\ are the SU(2)l and C/(l)y gauge coupling constants, and Yu,D,e,^ 
are the Yukawa matrices for up quarks, down quarks, charged leptons and neutrinos 



37 



respectively. The RGEs for Yd and the rest of the soft masses or trilinear terms are 
the same as in the MSSM. Below M, the RGEs are the same, except that the right 
neutrinos decouple, and hence Y u disappears from the equations. 

On the other hand, following ref.@ the amplitude for the process U — > Ijj can be 
written as 



T = e%m k ia a pqP{A L P L + A R P R )k 



At r, - 4 (C) 4- 4 (n) 



(75) 



A? 


e 1 


167r 2 m~ x 


Af 


e 1 


167T 2 m~ 


,(c,n) 
R 


— A L \L^R, 



M~- 

Cf AX C%xh{Ml-/ml x ) + C f AX C% x ^h{Mymi 



jh atL* j ti\j2 /„,2 \ i atL atR* 7 / i r- 



rrir 



(76) 

(77) 
(78) 



where M--(M^) is the chargino (neutralino) mass and m 2 x (m~ ) is the sneutrino 

(charged slepton) mass squared. 

The functions lip and Jx,2 are defined as follows: 

1 



h{r) 
hir) 

Mr) 
Mr) 



12(1 -r) 4 

1 ( 
2(1 -r) 31 

1 

12(1 -r) 4 
1 

2(1 - rf 



(2 + 3r — 6r +r +6rlnr), 
-3 + 4r - r 2 - 21nr), 
(1 — 6r + 3r 2 + 2r 3 — 6r 2 lnr) 
1 — r 2 + 2r lnr), 



Finally, 



= -92(O r ) a1 U x 



r L -a m 

u iAX — 92 



;{0 L ) A2 U X; 



(79) 
(80) 
(81) 
(82) 

(83) 
(84) 



\/2m w cos (5 

mii (0 N ) A3 U l Xii+3 }, (85) 
(0 N ) A3 U l Xti + 2(0 N ) A1 tan9 w U l Xii+3 }, (86) 



where A = 1, 2 and X = 1, 2, 3, and 

= -^%{[-(Oiv)A2 - {O N ) M tm0w)U l Xii + 

V 2 mvF cos p 

N l _ ^ 2 r m 'i 

m^cos/? 

where A = 1, ...,4 and X = 1,...,6. The different matrices in these equations are 
defined in Appendix B of ref . M . 
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